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1. Introduction 

The atmospheres of M stars are dominated by a small number of very strong 
molecular compounds (H2O, TiO, H2, CO, VO). Most of the hydrogen is 
locked in molecular H2, most of the carbon in CO; and H2O, TiO and VO 
opacities define a pseudo-continuum covering the entire flux distribution of 
these stars. The optical "continuum" is due to TiO vibrational bands which 
are often used as temperature indicators for these stars. These may be the 
depth of the bands relative to the troughs in between them; or the depth 
of the VO bands; or of the atomic lines relative to the local "continuum" ; 
or even the strength of the infrared water bands; all of these depend on the 
strength of the TiO bands and the amount of flux-redistribution to longer 
wavelengths exerted by them. Departures from LTE of the Ti I atom, and 
thus the concentration of the important TiO molecule, could, therefore, 
have severe and measurable consequences on the atmospheric structure and 
spectra of these stars. 

Due to the very low electron temperatures, the electron density is ex- 
tremely low in M stars; even lower than in novae and SNe. Collisions with 
particles other than electrons, e.g., the H2 or helium, are by far not as 
effective as electron collisions, both because of their smaller cross-sections 
and of their much smaller relative velocities. Therefore, collisional rates 
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Figure 1. Simplified Grotrian diagram of our Ti I model atom. All 395 levels and 5379 
primary (i.e., full NLTE) lines are shown but the 0.8 million secondary (approximate) 
NLTE lines have been omitted for clarity. 

which tend to restore LTE, could be very small in cool stars. This in turn 
could significantly increase the importance of NLTE effects in M stars when 
compared to, e.g., solar type stars. 

In this paper we discuss NLTE effects of Ti I in fully self-consistent 
models for a few representative M/Brown dwarf and M giant model atmo- 
spheres and spectra. 

2. Methods and Models 

In order to investigate if Ti I NLTE effects are important for the formation 
of cool star spectra, a full NLTE model calculation is required. This means 
that the multi-level NLTE rate equations must be solved for a number 
of species self-consistently simultaneously with the radiative transfer and 
energy equations, including the effects of line blanketing and of the molec- 
ular equation of state. For the purpose of this analysis we use the model 
code PHOENIX originally developed for the modeling of the expanding at- 
mospheres of novae and supcrnovac, and adapted to conditions prevailing 
in cool stars by Allard and Hauschildt (1995, hereafter AH95). PHOENIX 
(version 5.9) uses a spherical radiative transfer for giant models (loggi3.5), 
and an equation of state (EOS) recently upgraded from 98 to 195 molecules 
using the polynomial partition functions by Irwin (1998). For this work, the 
EOS as been extended to include: 

(a) TiO+ and ZrO+ (Gurvich & Glushko, 1982). 

(b) (and H^) to the EOS using the new partition function by Neale 
& Tennyson (1995). 

The opacity data base of Phoenix includes assorted f-f processes, H~ ab- 
sorption. Collision Induced Absorption (Borysow 1993), JOLA bands for 
CaH, VO, and FeH (AH95), dust absorption (experimental) and a direct 
line-by-line treatment for all available absorbers (inculding isotopes): 

(a) atomic lines from updated Kurucz CD#1 

(b) J0rgensen TiO & CN lines, H2O (opacity sampling tables) 

(c) Goorvitch CO lines 

(d) the inclusion of the IR roto- vibrational system of SiO using the detailed 
line list by Goorvitch. 

(e) hues by Neale & Tennyson (1995) 

(f) Kurucz CD#15 molecular lines 

(g) Miller & Tennyson H2O lines for high-res IR spectra 
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Figure 2. Ti I departure coefficients as function of optical dcptli for a NLTE model 
with Teff = 2700 K and log 5 = 5.0 (left hand side) and with with T^s = 4000 K and 
logs' = 5.0 (right hand side). 

(h) Improved water vapor lines by Viti et al (in preparation) as soon as 
they become available. 

For the strongest ca. 6 x 10^ atomic & molecular lines, we use detailed 
depth-dcpcndcnt Voigt profiles with improved damping constant compu- 
tation by Schweitzer et al (in preparation), and Gaussian profiles for an 
additional 10 x 10^ much weaker lines. In addition, we include ca. 2000 
photo- ionization cross sections (Vcrner & Yakovlcv 1995). Both the NLTE 
and LTE lines for atoms and molecules are treated with a direct line-by-line 
(LBL) method. This method is different from the classical opacity sampling 
(OS) approach in that we do not use pre-computed tables for the line opac- 
ity as a function of temperature using only fixed Gaussian line profiles. 
The LBL method dynamically selects the relevant lines from the master 
line lists at the beginning of each iteration and sums the contribution of 
every line within a search window to compute the total line opacity at ar- 
bitrary wavelength points. The latter is important in NLTE calculations in 
which the wavelength grid is both irregular and variable (from iteration 
to iteration due to changes in the physical conditions) . The LBL approach 
also allows detailed and depth dependent line profiles to be used during the 
iterations. However, in order to make the LBL method computationally fea- 
sible, modern numerical techniques, e.g., block algorithms with high data 
locality, and high-end workstations or supercomputers must be used. In 
the calculations we present in this paper, we have set the micro-turbulent 
velocity ^ to 2kms~^. We include LTE lines if they are stronger than a 
threshold V = xi/i^c = 10~^, where xi is the extinction coefficient of the 
line at the line center and Kc is the local b-f absorption coefficient. 

PHOENIX is a full multi-level NLTE code, i.e., the NLTE effects are 
included self-consistently in the model calculations. Hauschildt & Baron 
(1995) have extended the numerical method developed by Hauschildt (1993) 
for NLTE calculations with a very detailed model atom of Fe II. In this sec- 
tion we describe how we apply this technique to add the Ti I atoms to the 
list of NLTE species already available in the model calculations. 

2.1. THE TII MODEL ATOM 

To construct the Ti I model atom we have selected the first 34 terms of 
Ti I. We include all observed levels that have observed b-b transitions with 
log {gf) > —3.0 as NLTE levels where g is the statistical weight of the 
lower level and / is the oscillator strength of the transition. This leads 
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Figure 3. NLTE effects on the Ti ionization and molecule formation for the 
Tea = 4000 K, logs = 5.0 model. 

to a model atom with 395 levels and 5279 primary transitions treated in 
detailed NLTE. That is, we solve the complete b-f & b-b radiative transfer 
and rate equations for all these levels including all radiative rates of the 
primary lines. A Grotrian diagram of this model atom is shown in Fig. 1. 
In addition, we treat the opacity and emissivity for the remaining nearly 
0.8 million "secondary" b-b transitions in NLTE, if one level of a secondary 
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T^„=3200K, log g=1.0 Giant Model 



T.,=4000K, log g=1.0 Giant Model 
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Figure 4- NLTE effects on the Ti I line at Avac ~ 9641 A for two giant models. The model 
parameters are Tcs — 3200 K, \ogg — 1.0 for the left hand panel and Tes = 4000 K, 
log 5 = 1.0 for the right hand panel. The LTE spectrum uses the same model structure 
as the NLTE spectrum but with all departure coefficients set to unity. 

transition is included in the model. A detailed description of the method is 
given in Hauschildt Sz Baron (1995). 

Photo-ionization and collisional rates for Ti I are not yet available. Thus, 
we have taken the results of the Hartree Slater central field calculations of 
Reilman & Manson (1979) to scale the ground state photo-ionization rate 
and have then used a hydrogenic approximation for the energy variation 
of the cross-section. Although they are only very rough approximations, 
the exact values of the b-f cross-sections are not important for the opacities 
themselves which are dominated by known b-b transitions of Ti I and other 
species. They do, however, have an influence on the actual b-f rates but this 
remains unimportant for the computational method used in this work. 

While collisional rates are important in hotter stellar atmospheres with 
high electron densities, they remain nearly negligible when compared to the 
radiative rates for the low electron densities found in cool stars. We have 
approximated bound-free collisional rates using the semi-empirical formula 
of Drawin(1961). The bound-bound collisional rates are approximated by 
the semi-empirical formula of Allen (1973), while the Van Regemorter's 
formula (1962) was used for permitted transitions. 

A more accurate treatment of this model atom requires the availabil- 
ity of more accurate collisional and photo-ionization rates for Ti I. In the 
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T^„=2700K, log g=5.0 Dwarf Model 



T =4000K, log g=5.0 Dwarf Model 
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Figure 5. NLTE effects on the Ti I line at Avac « 964lA for two dwarf models. The 
model parameters are Tcb = 2700 K, log g = 5.0 for the left hand panel and Tefr = 4000 K, 
log 5 = 5.0 for the right hand panel. The LTE spectrum uses the same model structure 
as the NLTE spectrum but with all departure coefficients set to unity. 



present calculations we have neglected collisions with particles other than 
electrons because the cross-sections are basically unknown. Additional col- 
lisional processes would tend to restore LTE, therefore, the NLTE effects 
that we obtain in our calculations should be maximized. 



3. NLTE effects on Ti I line profiles in M Star Spectra 

We demonstrate the effects of NLTE on the formation of Ti I lines in 
Figs. 4 and 5. In general, the NLTE effects are smaller for the cooler models 
and they are larger for the giants than for the dwarfs at given effective 
temperature. In the optical spectral region, the changes caused by the Ti I 
NLTE line formation are very small and would be hard to observe due to 
the enormous crowding of lines in this spectral region. 

In the cooler models, the Ti I lines form deeper in the atmosphere in 
a region in which the radiation field is nearly Planckian and thus NLTE 
effects are very small. This is due to the enormous background opacity of 
TiO and other molecules. In the outer atmosphere of the cooler models the 
concentration of Ti I is much smaller than that of TiO, thus the effect of 
the large departures from LTE that we find in these regions on the Ti I line 
profiles is very small. 
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In the hotter models, NLTE effects, in particular on the near-IR lines, 
are much larger. In these models the line forming region of the Ti I lines 
is inside the region in which the departures from LTE arc significant. In 
addition, the TiO opacity is relatively smaller than in the cooler models. 
For the Ti I line at Avac ~ 964lA, NLTE effects make the core of the line 
deeper than the LTE model predicts. This is due to line scattering which 
removes photons from the line core and re-distributes them into the line 
wings. The same effects are present in most Ti I lines, but there are a few 
exceptions for which the line is weaker in NLTE than in LTE. Abundance 
determinations of Ti, or likely all metals, from near-IR or IR lines should 
therefore include NLTE effects wherever possible. 

For the model parameters that we have considered so far, the effects of 
Ti I NLTE on the TiO bands is very small. This seems to be due to the fact 
that the line forming region of TiO (around Tstd ~ lO""' is dominating the 
line forming region of the Ti I lines (due to the larger opacity of the TiO 
molecule compared to Ti I). Therefore, the TiO lines form in a region in 
which the Ti I atom is basically in LTE and the Ti I NLTE effects that are 
important at smaller optical depths cannot affect the TiO lines significantly. 
The situation is different for the TiO^ molecule which forms from Ti II and 
O I (cf. Fig. 3). This molecule is very sensitive to the Ti I NLTE effects 
and its lines would be helpful indicators for Ti I NLTE effects. 

In further investigations it will be very interesting to look for NLTE 
effects in the TiO line formation itself. This is certainly feasible with modern 
numerical techniques once adequate data for TiO are available. 
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